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ABSTRACT 

A naive extension to MOND of the standard computation of the two-body relaxation 
time t2h implies that t2b is comparable to the crossing time regardless of the number 
N of stars in the system. This computation is questionable in view of the non-linearity 
of MOND's field equation. A non-standard approach to the calculation of t2b is de- 
veloped that can be extended to MOND whenever discreteness noise generates force 
fluctuations that are small compared to the mean-field force. It is shown that this 
approach yields standard Newtonian results for systems in which the mean density 
profile is either plane-parallel or spherical. In the plane-parallel case we find that in 
the deep-MOND regime <2b scales with N as in the Newtonian case, but is shorter 
by the square of the factor by which MOND enhances the gravitational force over its 
Newtonian value for the same system. Near the centre of a spherical system that is in 
the deep-MOND regime, we show that the fluctuating component of the gravitational 
force is never small compared to the mean-field force; this conclusion surprisingly 
even applies to systems with a density cusp that keeps the mean-field force constant 
to arbitrarily small radius, and suggests that a cuspy centre can never be in the deep 
MOND regime. Application of these results to dwarf galaxies and groups and clusters 
of galaxies reveals that in MOND luminosity segregation should be far advanced in 
groups and clusters of galaxies, two body relaxation should have substantially modified 
the density profiles of galaxy groups, while objects with masses in excess of ^ IQMq 
should have spiralled to the centres of dwarf galaxies. 

Key words: gravitation - galaxies: dwarf - galaxies: haloes - galaxies: kinematics 
and dynamics 



1 INTRODUCTION 

Milgrom (1983) proposed that the failure of galactic rotation 
curves to decline in Keplerian fashion outside the galaxies' 
luminous body arises not because galaxies are embedded 
in massive dark haloes, but because Newton's law of grav- 
ity has to be modified for fields that generate accelerations 
smaller than some value ao. Bekenstein & Milgrom (1984; 
hereafter BM84) proposed the non-relativistic field equation, 
eq. (|5J below, for the gravitational potential $ that gener- 
ates an appropriately modified gravitational acceleration g 
of a test particle through 

g^-V*. (1) 

A considerable body of observational data now supports 
this theory of modified Newtonian gravity (MOND) - see 
Sanders & McGaugh (2002) for a review. The big problem 
with MOND is our inability to derive the MOND field equa- 
tion as the low-energy limit of a Lorentz covariant theory. 
This inability is unfortunate in two respects. First it makes 



it impossible to determine MOND's predictions for gravita- 
tional lensing experiments, or any observation that involves 
relativistic cosmology, particularly observations of the CMB. 
Since these are the areas in which the competing Cold Dark 
Matter (CDM) theory has been most successful, the lack of 
a Lorentz covariant form of MOND makes a fair confronta- 
tion between MOND and CDM impossible. Another reason 
to regret this lack is that there are tantalizing hints that the 
characteristic acceleration ao ~ 2 x 10~* cm s~^ that lies 
at the heart of MOND is connected to the requirement for 
a non-zero cosmological constant A in Einstein's equations: 
ao and A ~ 3(ao/c)^ may be two aspects of a single physical 
process associated with the unknown small-scale structure 
of space-time (Milgrom 2002). 

Clearly it is worthwhile to pursue these ideas regard- 
ing a putative Lorentz covariant form of MOND only if 
the low-energy theory that we already have accounts for 
all available data. In this paper we show that MOND pre- 
dicts two-body relaxation times for systems in which g <^ ao 
that are shorter than those given by Newtonian theory by a 
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factor ~ (g/ao)^. We show further that in MOND dynam- 
ical friction operates on a timescale that is shorter than in 
Newtonian dynamics with dark matter by a factor ~ g/ao- 
The shortness of the dynamical friction timescale has obser- 
vationally testable predictions for the dynamics of systems 
such as dwarf galaxies, groups and clusters of galaxies that 
are in the deep MOND regime. 



2 ORDER-OF-MAGNITUDE ANALYSIS 

BM84 replace Poisson's equation for the potential in terms 
of the density p by 



V- [At(lV'l'l/ao)V$] =A-kGp, 



(2) 



to be solved subject to the boundary condition |V0| 
for |x| — > oo. The function p,(x) is required to have the 
behaviour 



X for X ^1, 
for X 3> 1, 



(3) 



but the detailed manner in which p moves between these 
limits is currently constrained by neither observational data 
nor theory. 

We shall be concerned with the 'deep MOND regime' 
in which ii{x) ~ x. In this limit equation ^ can be written 



|V$|V$ 



do 



0, 



(4) 



where «1?n is the Newtonian potential generated by the given 
density distribution. This equation implies that the differ- 
ence between the two terms in the square brackets is equal to 
the curl of some vector field. BM84 show that when the den- 
sity distribution is spherical, planar or cylindrical, the curl 
vanishes. It then follows that the acceleration go in MOND 
is related to the Newtonian acceleration gN by 



(5) 



We now present a heuristic derivation of the two-body 
relaxation time of a system that is in this regime. Our analy- 
sis is modelled on the standard Newtonian derivation of the 
two-body time (e.g., §4.1 of Binney & Tremaine 1987, here- 
after BT). This derivation is straightforward to follow and 
illuminates the basic physical principle that causes two-body 
relaxation to be fast in MOND. It is open to the criticism, 
however, that it ignores the inherent non-linearity of the ba- 
sic equation Q. Consequently, in the following section we 
rederive the Newtonian relaxation time by a very different 
technique that carries effortlessly over to the case of MOND. 

Milgrom (1986; hereafter MSG) shows that in the deep 
MOND regime equation Q causes a force F to act between 
two isolated point masses mi and m2, which can be written 



Gm\m2 r / rn2 r 

^ = 5 / ' — 

\ mi ro 

where 



G{mi + m2) 
do 



: 8.1 10' 



16 mi+m2 



(6) 



(7) 



The function / in equation Q can be calculated numerically 
as a function of r/ro for a given mass ratio. In particular, 
numerical/asymptotic solution of equation (|5J shows that: 



• for r/ro < 1, / ~ 1 independently of the mass ratio; 

• for 1712 ^ mi, the acceleration of m2 is given by equa- 
tion Q, while the acceleration of mi is obtained from the 
conservation of linear momentum migi + m2g2 = 0; 

• finally, for nearly equal masses in the deep-MOND 
regime, 



mim2 



Gm\ra2 
rro 



(8) 



\/mr+m^ r 

According to the numerical calculations of M86, this for- 
mula does not err by more than 20% for any mass ratio. 
Note that its asymptotic behavior is correct when one of 
the masses is vanishingly small or much more massive than 
the other. 

For ordinary stars, ro is much less than the mean inter- 
stellar distance in galaxies, so nearly all stellar interactions 
are in the deep-MOND regime if the local gravitational mean 
field g of the whole galaxy is. If, by contrast, g>ao, then all 
interactions, even ones that are individually weak, conform 
to Newtonian theory - this result is sometimes called the 
"external field effect" (BM84). An immediate consequence 
of this result is that MOND predicts standard two-body re- 
laxation times for objects in which the mean field exceeds 
ao, even though individual stellar interactions generate ac- 
celerations much smaller than ao. 

Consider then a stellar encounter in a system that is in 
the deep-MOND regime. For simplicity we assume that both 
stars have the same mass m. Let the encounter be charac- 
terized by impact parameter b> ro and asymptotic relative 
velocity V. Using equation ^ and the impulse approxima- 
tion, we conclude that the encounter changes the velocity of 
each star (in the direction perpendicular to V) by an amount 



2b F{b) 
V ~ 



2Gm 



(9) 



Remarkably, this formula does not contain the impact pa- 
rameter b: in the deep-MOND regime, the deflection is inde- 
pendent of the impact parameter. The steady accumulation 
of such velocity changes causes the star's velocity to exe- 
cute a random walk. If the system contains stars and has 
half-mass radius R, per crossing time the star experiences 

5n = ^2Tvbdb (10) 

encounters with impact parameter in {b -f db,b). Adding in 
quadrature the velocity changes from successive encounters 
on the assumption that they are uncorrelated yields a cumu- 
lative change in the square of the stellar speed per crossing 
time 



(At;i)c 



N f2Gm\ 



R 



Abb = 2N(^\ 
\roV J 



(11) 



i?2 V roV J 

The ratio of the two-body time to the crossing time is thus 



t2h 



An 



tc 



(A^5 



2G^AP 



(12) 



where M — Nm is the system's mass and vtyp a typical 
stellar velocity. Finally we use equation Q and the funda- 
mental MOND equation vfyp = GMao to eliminate ro and 
Vtyp from this equation, and find 

t2h = Across. (13) 
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Thus this naive calculation implies that in the deep-MOND 
regime the two-body relaxation time is comparable to the 
crossing time regardless of the number of stars A*'. The phys- 
ical origin of this result is clear: equation Q states that in 
the deep-MOND regime distant encounters produce much 
larger deflections than in the Newtonian case, and even in 
the latter case distant encounters make a large contribution 
to the relaxation rate. 



3 RELAXATION IN A UNIFORM FIELD 

The derivation of t2h that we gave in the last section is 
objectionable on two grounds. First, in MOND all two-body 
orbits are bound because the potential is asymptotically log- 
arithmic. So our use of the impulse approximation is highly 
suspect. Second, the derivation assumes that the effects of 
encounters can simply be added. The dominant encounters 
are those at impact parameters comparable to the half-mass 
radius. Each such encounter lasts of order a crossing time, so 
the many encounters that contribute to (Ai;^)cross in equa- 
tion mil l occur simultaneously. In the deep-MOND regime 
the field equation ^ is highly non-linear, and it is far from 
clear that the effects of different encounters can be simply 
added. In this section we derive t2b by a different approach 
that is not open to this objection. 

The underlying physical idea is that two-body relax- 
ation is driven by fluctuations in the gravitational potential 
of a system that is in virial equilibrium. The fluctuations 
are generated by Poisson noise. We obtain t2h by decompos- 
ing the fluctuations into different spatial frequencies, and 
summing over frequencies rather than over encounters. We 
shall find that with this approach, the calculation of t2h for 
MOND differs very little from the corresponding Newtonian 
calculation. Hence we now rederive the Newtonian relax- 
ation rate with the new formalism, to demonstrate that it 
produces the familiar result, and to pave the way for the 
calculation of t2h in MOND. 



3.1 Newtonian relaxation 

We wish to consider the case in which the underlying system 
generates a uniform (Newtonian) gravitational field g. Such 
a field is generated by an infinite sheet of constant density. 
So we consider the effect that density fluctuations in this 
sheet have on a star that is located distance z from the 
sheet. Let x be a two-dimensional vector of coordinates in 
the plane of the sheet. Then at z 7^ Laplace's equation has 
solutions (e.g. BT §5.3.1) 



$(x) 



d^k$(k)e-*l"'exp(ik.x) 



(14) 



where $ is an arbitrary function, k — |k|, and d k = 
d^k/(27r)^. If the surface density of the sheet is 

E(x) = y d^kE(k)exp(ik.x), (15) 

then an application of Gauss's theorem shows that 
,S(k) 



Differentiating 11411 with respect to x and integrating with 
respect to time, we calculate a component of velocity parallel 
to the sheet that the fluctuating density induces in our test 
particle: 



Vx{r) = - / dt — 
Jo 9x 

= 27rGy dt J d^k^S(k,t)e-'=l"'exp(ik.x). (17) 
Squaring Vx and taking an ensemble average, we have 

X / dt / dt' exp[i(k-|-k').x] /E(k,t)E(k',t')V(18) 



Let Xq, denote the location within the sheet of star a. Then 
E(x,t) =m^(5[x<,(t) -x] (19) 



E(k,t) ^ exp[— ik.Xa(t)]. 



(20) 



For the moment we assume that the velocity of star a 
lies within the sheet. Then to a sufficient approximation 
Xa(t') = yia{t) + (t' — t)vQ and we have for uncorrelated 
stars 



E(k,t)E(k',t')) = 



X ^^exp[-ix^.(k-f k')]exp[i(t'-t)k'.vj^ . (21) 

Since the velocities of stars are independent of their posi- 
tions, the ensemble average above can be expressed as the 
product of two ensemble averages, one over x^ and the other 
over Vo, . With K = k -I- k' we have 

(exp(-ix<,.K)) = J ^exp(-ix.K) = i^5(K), (22) 

where A is any large area of the sheet, so 

^ {exp(-iXc.K)) = {2nfn5(K), (23) 

a 

where n is the number of stars per unit area. For the other 
ensemble average we have for a Maxwellian distribution of 
velocities 



(exp[i(t' - t)k'.VQ]) = 



27ra" 



e " exp[i(t' — t)k'.v] 



(24) 



= exp|^— — t) fc (T 

Substituting equations 12111 . I23II and 12411 into equation list 
and integrating over t' and t, we find for r 2> fef 



d^ke-""i"i^ 



(27r)^/"G^m"n 
4\z\a 



(25) 



<I>(k) = -2-kG- 



k 



(16) 



From equation l|25|l we can recover the standard expres- 
sion for the Newtonian diffusion coefficient by summing the 
contributions from many sheets. If p is the (homogeneous) 
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mass density due to stars, then pdz = nm, and the overall 
diffusion coefficient is 



2cr 



z 



pm 



2a 



In A, (26) 



(27) 



with A = imax/^min- This diffusiou coefficient may be com- 
pared with half the value of D{Avj_) in equation (8-68) of 
BT. Taking the limit X ^ of small test-particle velocities 
we find 

2(^g(T))^ ^ 6n 
tD{AvI) ^ 8 ■ 

Thus this derivation agrees with the classical one to within 
the uncertainties inherent in either approach. The weakest 
part of the present derivation is the assumption that the 
velocities of field stars lie within planes z — constant. Re- 
laxing this assumption would reduce the auto-correlation of 
each sheet's surface density below that given by equation 
12H and introduce correlations between the densities of dif- 
ferent sheets. When the change in z during an encounter 
is small compared to the distance of the particle from the 
point of observation, little will have changed physically from 
the case of constant z, so the new correlations will almost 
exactly compensate for the lost contribution to the autocor- 
relation. This argument suggests that the error introduced 
by confining particles to planes of constant z is not large, as 
is also indicated by the agreement between our value of the 
diffusion coefficient and that obtained in the standard way. 



3.2 Relaxation in MOND 

So long as the total number of stars in the system is large, 
and we avoid special points of symmetry such as the sys- 
tem's centre, the ffuctuations in the gravitational field are 
small compared to the MOND mean field g. Hence we may 
linearize the field equation ^ around g and solve a lin- 
ear field equation for the component of the potential that 
drives relaxation. As in the Newtonian case, we decompose 
the perturbing density field into sheets normal to the mean 
gravitational field. We find the potential ffuctuations due to 
each sheet, and add the effects of the sheets. 

M86 shows that when the field equation JJJ is perturbed 
about a uniform gravitational field go, the first-order per- 
turbations to the potential and density are connected by 



Oo 

90 



4-kGp, 



(28) 



where V is the full three-dimensional Laplacian operator. 
Thus the solutions 11411 to Laplace's equation must be re- 
placed by 



$(x) 



d^k$(lc)e 



exp(ik.x). 



and equation (I16I I becomes 
,S(k) 



<E.(k) = - 



2tvG- 



k 



(29) 



(30) 



The Newtonian calculation carries over with two substitu- 
tions: \z\ — ^ \z\l^2 and m — > [ao/ ^2gQ)m. Equation I25II 
for the contribution of a single sheet to the diffusion coeffi- 
cient becomes 



{vI{t)) _ {2-Kfi^G^n?n ag 



(31) 



4|z|a V2flo" 
so the diffusion coefficient is larger than in the Newtonian 
case by a factor aQ/(^2t;o). 

Consider now two systems that contain identical distri- 
butions of stars in phase space, but differ in the way gravity 
works: in one system gravity is Newtonian, while in the other 
it is described by MOND. In the Newtonian case we augment 
the gravitational field of the stars with a fixed background 
field to ensure overall dynamical equilibrium. Then the two- 
body relaxation times in the two systems are in the inverse 
ratio of their diffusion coefficients, which from equations 1251 1 
and lEB 

is 



,M 

''2b 



(32) 



This result will remain true when the fixed background field 
in the Newtonian system is replaced by the field generated 
by a distribution of DM particles provided individual DM 
particles are much lighter than stars - the particles then 
make negligible contributions to the fiuctuations in the over- 
all gravitational acceleration^ . 

In the deep-MOND regime, equation Q holds, so elim- 
inating ao from equation 1321 . the ratio of two-body times 
becomes 



''2b 

where 



9l 



V2 



(l + 7^)2 



7^; 



(33) 



(34) 



is the ratio of the DM to stellar mass in the Newtonian 
system. 

3.3 Dynamical friction 

Equation 13311 has important implications for the magnitude 
of dynamical friction in MOND because the dynamical fric- 
tion experienced by a body as it moves through a popula- 
tion of background objects that are in thermal equilibrium, 
is proportional to its diffusion coefficient in velocity space. 
This relation can be verified for coefficients calculated un- 
der the assumption of standard Newtonian gravity, but, as 
Chandrasekhar (1943) pointed out, it follows from the gen- 
eral principles of statistical mechanics. This fact ensures that 
the relation is valid regardless of the law of gravity. 

In the local approximation, the Fokker-Planck equation 
for the evolution of the distribution function / of a popula- 
tion of 'test' objects can be written 

At av ■ ' 

where the flux S of stars in velocity space is [BT eq. (8-57)] 



(35) 



1 



d 



2 ^ d% 



(/A,). 



(36) 



In a Ne-wtonian system -with field particles of two different 
species, with masses m* and moM and densities p» and poMi 
the two body relaxation time of a test particle of mass ma is 
*2b ^ l/(m.P« +mDMPDM)- 
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Here Di = D{Avi) and Dij = D{Avi Avj) are the usual 
first- and second-order difi'usion coefiicients. Let the scat- 
tering objects have mass mf and be in thermal equilib- 
rium at inverse temperature f3 — {mfa^)~^ . Then the prin- 
ciple of detailed balance implies that S will vanish when 
/ ~ exp{—l3H), where H = ma(|i'^ + $) is the Hamiltonian 
that governs the motion of the test objects. In a frame in 
which Dij is diagonal it then follows that 



Di = —^(3ms,ViDii + 



1 A 



/ ma 2 d 
\2mi 



2 a ^ \ 



(37) 



(no sum over repeated indices). For Vi < a the term with 
the logarithmic derivative is of order unity - in the case of a 
Maxwellian distribution of scatterers it evaluates to —3/10 
at Wi = [BT, eq. (8-65)]. Hence for ma Z 4mf equation 13711 
implies that Vi decays exponentially in a characteristic time 
tfric that is related to the two-body relaxation time by 

tfric _ 2mf ^ 

t2b ma ' 



(38) 



this result is consistent with the Spitzer (1987) relation 
ifric/t2b ~ 2mf/(mf -I- ma), when ma^J 4mf. From this re- 
sult and eq. 1331 it follows that in MOND the friction time 
is reduced by a factor (1 -f 71)^/^2 over the value it would 
have in a Newtonian system with the same stellar mass and 
a fixed auxiliary gravitational field. 

Whereas the diffusion coefficients were unchanged when 
the fixed gravitational field was replaced by the field of 
swarms of low-mass DM particles, this replacement enhances 
dynamical friction by a factor (1 -f TV). Hence^. 



2b 

and 

,M 

'•fric 

j-N 

'•fric 



1 + n' 



(39) 



(40) 



4 EXTENSION TO SPHERICAL SYSTEMS 

Given that stellar systems are frequently approximately 
spherical, and rarely have the plane-parallel symmetry that 
we assumed in the last section, we try to adapt the preced- 
ing calculation to a spherical system. The appendix derives 
an expression for the Newtonian diffusion coefficient expe- 
rienced by a star that is stationary near the centre of a 
spherical system. 

If we are to use perturbation theory to carry this New- 
tonian analysis over to the deep-MOND regime, the unper- 
turbed acceleration go should be non- vanishing and less than 
oo at the location of the test star. For simplicity the sys- 
tem's density profile is taken to be a power law in radius 

^ A system with comparable masses of luminous and DM vio- 
lates equation 1381 because the stars and DM particles are not in 
thermal equilibrium with one another. The standard treatment 
of dynamical friction shows that for a test particle of mass ma 
t^j^ oc l/[(ma -I- m«)p« -I- (ma -I- mDM)PDM], and from Footnote 
1, when ma Z of m* and mDM, and m« > m^M equation (35) is 
reobtained. 



p r"'', and the requirement that at all radii go/ao be 
of order but less than unity then implies 7 = 1. We show 
that in this case the acceleration caused by the fluctuations 
becomes comparable to go as one approaches the centre, 
and the linearized field equation ceases to be valid. We infer 
from this result that in the deep-MOND regime the centres 
of all systems must be homogeneous, since non-linear fiuc- 
tuations in the acceleration will soon disrupt the cusp in the 
density profile that generates the assumed constant acceler- 
ation go- Unfortunately, perturbation theory cannot be used 
to calculate the central relaxation time of a system with a 
constant-density core. 

Thus this analysis does not lead to a value of t2h for 
systems in the deep-MOND regime, but it does suggest that 
at the centres of these systems two-body relaxation is very 
much more rapid than in the Newtonian case. 



5 ASTROPHYSICAL APPLICATIONS 

We now apply the results of Section 3 to stellar systems 
that in Newtonian dynamics would turn out to be DM dom- 
inated, namely to dwarf galeixies, to galaxy clusters and 
to galaxy groups: all these systems are in the deep-MOND 
regime. 

Equation (8-71) of BT gives the Newtonian two-body 
relaxation time for a system with a given velocity dispersion 
a and density of scatterers p«. As in BT, we substitute into 
this equation values of a and p* appropriate to the system's 
half-mass radius. We obtain a result for the reference relax- 
ation time that is analogous to equation (8-72) of BT but 
different because now ~ 0.4G'M»(1 + TZ)/rh on account 
of the presence of DM. With A = 0.4iV(l + 7^) we have 



t' 



0.66 Gyr 



M, 



In A m \pcj I^IO^Mq 



1/2 



(41) 



while equation 1331 now implies that the relaxation time for 
MOND is 



t.-t, ~ 0.9 Gyr — — 

In A m \ pc 



3/2 



M, 



105 Af0 



1/2 



(42) 



5.1 Dwarf galaxies 

For a dwarf galaxy such as Draco, AI, ~ 2.610^Mq, rh — 
200/i~^ pc, and 7^ ~ 100 (see. e.g., Mateo 1998, Kleyna et al. 
2002), so the reduction factor is enormous (~ 10*). However, 
tjlji — 10^ Gyr, so tj^ is still is slightly longer than the Hubble 
time. 

Since in MOND the dynamical-friction time for inspi- 
ralling of an object of mass ma is shorter than the two-body 
time by ~ Mo/ma, any object significantly more massive 
than a star will have spiralled to the centre of a dwarf galaxy. 
In particular, any black holes with masses J; IOMq should 
have collected at the centre. Since globular clusters have 
masses > lO^M©, they are liable to spiral to the galaxy cen- 
tre even in Newtonian dynamics. In MOND they spiral in 
on essentially a dynamical time. Consequently, the posses- 
sion of globular clusters by a dwarf galaxy that is in the 
deep-MOND regime would be problematic for MOND. In- 
terestingly, the dwarfs with large values of TZ, namely Draco, 
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UMi, Carina, and Sextans, have no globular clusters listed 
in Table 8 of Mateo (1998). 

5.2 Galaxy clusters 

If Mt is the total cluster mass, we have M« ~ 0.05AfT, 
Mgas ~ O.ISMt, Mdm ^ 0.8Mt (see, e.g., Bohringer 1996), 
so 7^ ~ 4. Adopting N ~ 100, ~ 1.7 10^''' M©, rh ~ 
580 kpc, and a = 1000kms"S we find ~ 48Gyr and 
C~2.7Gyr. 

In MOND dynamical friction will cause a galaxy of 
mass nia to spiral to the cluster centre on a timescale 
~ 5{rn/ma) Gyr, where m is the mass- weighted mean galac- 
tic mass. It follows that in MOND luminosity segregation 
should be well advanced within galaxy clusters. 

5.3 Galaxy groups 

For a typical galaxy group, ~ 5, M, ~ 5.2 lO^Af©, 7^ = 
16, a ~ 200kms"\ and rh ~ 380 kpc. Hence - 70 Gyr 
and ~ 340 Myr. If one (improperly!) applied the concepts 
of Newtonian stellar dynamics, one would conclude that in 
MOND galaxy groups should have already evaporated be- 
cause in the Newtonian case the evaporation time is ~ lOtrh- 
The logarithmic asymptotic form of MOND's two-body po- 
tential implies that it is impossible to escape from a perfectly 
isolated MOND system''. However, the shortness of for 
galax;y groups suggests that relaxation will have significantly 
increased the central densities of groups. Luminosity segre- 
gation will proceed on a dynamical timescale. 



6 DISCUSSION AND CONCLUSION 

A naive extension to MOND of the standard derivation of 
the two-body relaxation time t^h leads to the conclusion that 
t2b is comparable to the crossing time for any value of the 
number A^" of stars in the system. The derivation is open 
to objection on at least two grounds, so we have developed 
an approach to the calculation of the Newtonian relaxation 
time that can be straightforwardly adapted to MOND. This 
focuses ou the effect of density fluctuations, which must be 
small in the limit of large A'^, and calculates the correspond- 
ing potential fluctuations by linearizing MOND's field equa- 
tion around a uniform background field go ^ ao. This anal- 
ysis reveals that in MOND t2b scales with in the same 
way as it does in the Newtonian case, but is smaller by a 
factor ~ (1 -|- 71)^, where TZ is the ratio of the apparent DM 
and stellar masses. 

We show that the timescale tfric on which dynamical 
friction causes an object of mass rria to spiral in through a 
population of scatterers of mass rrif is t^?;^ = (2mf/ma)t2b 
in the case of MOND, while the Newtonian time is tf^j^ ~ 
{l+TZ)t^i^, provided individual DM particles have negligible 
masses. 

^ We note that it is generally supposed that the galaxy disks 
are made of shattered star clusters and associations: while tidal 
destruction will still work, the failure of evaporation in outlying 
parts of galaxies might lead to the survival of surprising numbers 
of low-mass star clusters. 



Application of these results to DM-dominated systems 
shows that at the half-mass radius, the two-body relaxation 
times for MOND of both dwarf galaxies and clusters of 
galaxies are of order the Hubble time. Consequently, MOND 
predicts that objects in these systems that have masses more 
than ~ 10 times the mean mass will have spiralled to the sys- 
tems' centres. In particular, any black holes formed in core- 
collapse supernovae should have collected at the centres of 
dwarf galaxies, and luminosity segregation should be far ad- 
vanced in clusters of galaxies. The globular clusters of dwarf 
galaxies should spiral inwards on a dynamical timescale, so 
the existence of globular clusters in DM-dominated dwarf 
galaxies would be problematic for MOND. 

In MOND the two-body and dynamical-friction 
timescales of groups of galaxies are only ~ 300 Myr. A New- 
tonian system with such a short relaxation time would have 
evaporated completely, but evaporation is probably impos- 
sible in MOND. However, in MOND relaxation will surely 
have significantly modified the density profiles of galaxy 
groups, and induced substantial luminosity segregation. Fur- 
ther work is required to determine whether such evolution 
is compatible with observation. 

The modification of Poisson equation given by equation 
(2) is only one possible formulation fo MOND. We have not 
considered relaxation in modified inertia theories, which give 
very similar predictions for rotation curves and global mass 
discrepancies (Milgrom 2002). These theories could yield 
very different predictions for relaxation phenomena. 
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APPENDIX A: RELAXATION AT THE 
CENTRE OF A SPHERICAL SYSTEM 

Al Newtonian case 

We consider the stochastic acceleration of a particle that is 
initially stationary at the centre of a spherical system. We 
expand the system's gravitational potential, generated by 
the (random) density field 
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pi 



= ^ ma(5(x - x„) 



(Al) 



in spherical harmonics. At the origin, only the dipole term 
makes a non-vanishing contribution to the force on a parti- 
cle. Moreover, near the origin the force contributed by the 
Ith multipole varies as r'~^ times a factor that changes sign 
each time the particle passes the origin for even /. Hence av- 
eraged along a trajectory through the orgin, the force from 
multipoles with / > 1 is smaller by that from the dipole term 
by a factor that vanishes with the square of the apocentric 
radius. We therefore concentrate on the dipole term, which 
is (BT §2.4) 



where 



(t>)r / Aap-im{a) 



(A2) 



pim(a)= / d?Q.Y^* p{a,6,(j)) 



= -!--|- / d^x Yi™* mQ5(x - Xa) 
da J ^ — ' 



(A3) 



where the sum over a is extended to particles inside the 
sphere of radius a. Differentiating IIA21 with respect a 
generic direction (which without loss of generality we as- 
sume to be z), and evaluating the obtained expression at 
r = 0, we have 



9$ 



Vz = ^ l'^*^ / dapio(a) 



(A4) 



because the contributions from m — ±1 vanish. Integration 
with respect to time gives 

Vz{r) = fyrG dt dapio(a). (A5) 
Jo ^0 

Squaring and taking the ensemble average we find 
(^') = JJ dada'^ dtdt' (pio(a,t)pio(a',t')) -(AG) 

The expectation value in the integrand is non-negligible only 
for |a — a'l and |t — sufficiently small. With IIA3t we have 
in the case that all particles have equal masses 



C{a,t,t') = J da' (pio(a, t)pio(a',t')) 
d f 1 d / 



a'^ da I a'^ da' ^—^ 

a, 13 



(Y;\a.t)yi'w,t')) 



a 

where we have assumed that the particles are mutually un- 
correlated. If a is the characteristic tangential velocity of 
particle a, the expectation value in this equation vanishes 
after a time of order raja- Therefore we replace the inte- 
gral over t' of C{a,t,t') by a/a times half its peak value. 



To estimate the latter we observe that by the Monte-Carlo 
theorem 



1 = / d^fl\Y"f 



Hence 

dt' C{a,t,t' 



N 

^2^in(a)i . 



da ( (T ) 



Svra^ da 

Inserting this expression into ijAB^ we have 

(vzi-r)) _ 2n{Gmf [ d{N/a) 



D 



9cr 



(A8) 



(A9) 



(AlO) 



If the stellar density p/m is a constant from some smallest 
radius out to some maximum radius rmax and then zero, then 
A*' — ^■Ka?{p/m) and the diffusion coefficient D becomes 



^ 167r G mp , 

D — In 

27a 



(All) 



This should be compared with half the value of 
limx-.oI>(Aui) in eq. (8-68) of BT. One finds 



2D 



= 1.75 



(A12) 
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As in the case of plane-parallel geometry we obtain a diffu- 
sion coefficient that is larger by a factor < 2 than its con- 
ventional value. 

A2 Case of the deep-MOND regime 

For the reasons given in the main text we assume that the 
system's density profile obeys the power law p ~ 1/r, which 
yields a constant radial acceleration go — ^JooGMo/iq Cr, 
where Mo is the mass interior to the fiducial radius ro. 

The perturbed potential $ satisfies eq. (7) of M86 with 
Ho = go /do and Lo = 1 because we are in the deep-MOND 
regime. Hence 

\^(r^'^)-^^^il,.Gp (A13) 
or \ or J r^ ao 

where C is the angular part of the Laplacian. We seek the 
dipole term in the expansion of <1> in spherical harmonics 



(A14) 



i = m=-i 

which satisfies 

— ^V^—] 5- = — 27rGpio. (A15) 

or \ or J ao 

The Green's function u(r, r) for this equation satisfies 

1 d ( od' 



dr 



r — = —2nGd{r, r ) 

\ or J ao 



(A16) 



From the power-law solutions to the homogeneous equation 
it follows that 



/^ = / A{r')r"-'/^ for r < r' 
^ ' ' \ B(r')r-'^-i/2 forr>r'. 



(A17) 



where A{r) and B{r) are functions to be determined from 
the rhs of eq. IIA16II and a = \/5/2. Consequently, near 
the origin the dipole term varies as r^^^~^'>/^ ^ ,.0-62 ^^-^^ 
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the acceleration to which it gives rise diverges as r~ . 
Hence no matter how small the fluctuations in p are, the 
linearization of ((2J breaks down sufficiently near the origin. 
We conclude that sufficiently near the origin the fluctuating 
part of the field is comparable to, or larger than, go - In these 
circumstances the fluctuations will soon disrupt the cusp in 
the density distribution that generates go. It seems therefore 
that in the deep-MOND regime the centres of all systems 
must be homogeneous. Unfortunately, perturbation theory 
cannot be used to calculate the central relaxation time of 
such a system. 



